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Abstract 

This paper considers the gain-scheduled leader-follower tracking control problem for a parameter 
varying complex interconnected system with directed communication topology and uncertain norm- 
bounded coupling between the agents. A gain-scheduled consensus-type control protocol is proposed 
and a sufficient condition is obtained which guarantees a suboptimal bound on the system tracking per¬ 
formance under this protocol. An interpolation technique is used to obtain a protocol schedule which is 
continuous in the scheduling parameter. The effectiveness of the proposed method is demonstrated using 
a simulation example. 
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1 Introduction 

In recent years, the topic of cooperative control has attracted much attention. The objective of the 
cooperative control problem is to propose distributed control laws to achieve a desired system behavior [1]. 
A related problem is that of synchronization of complex dynamical systems where all the components are 
controlled to exhibit a similar behavior by interconnecting them into a network [2, 3]. 

There are several approaches to the synchronization problem for complex systems consisting of many 
dynamic subsystems-agents. In the average consensus problem, which has received considerable atten¬ 
tion in the last decade [4, 5], the objective is to synchronize all the agents to a common state. Another 
approach is to employ a suitable internal model for the system. For instance, it is shown in [6] that the 
existence of such an implicit internal model is necessary and sufficient for synchronization of a system 
of heterogeneous linear agents considered in that paper. However, the internal model approach does not 
directly address the overall system performance, since each system is controlled to follow its own internal 
model dynamics. In general, it may be difficult to determine an internal model that guarantees a good 
synchronization performance. Yet another approach is to designate one of the agents to serve as a leader, 
and design interconnections within the system so that the rest of the system follows the leader [7]; also, see 
[8] for a recent example. This idea leads to the leader-follower tracking problem, which is the main focus 
of this paper. 

A common feature of many papers that consider the leader-follower problem is that the dynamics 
of agents are usually assumed to be dynamically decoupled [22, 23, 24, 25, 26, 27]. In many complex 
systems, however, interactions between subsystems are inevitable and must be taken into account [11]. 
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Examples of systems with dynamical interactions between subsystems include spacecraft control systems 
and power systems [9]. For example, it was noted in [10] that weak uncertain couplings between generators 
is one of the reasons for dynamic instability in power systems and, therefore, they should not be neglected. 
This motivates us to consider the leader-follower tracking problem for interconnected systems. While we 
do not consider a specific application, our approach makes a step towards using consensus feedback for 
synchronization of such systems, compared to other contributions in the area of networked control systems 
which do not address the presence of interconnections. 

In this paper, we are concerned with the leader-follower tracking problem for interconnected systems 
that depend on a time varying parameter. Analysis and control of linear parameter varying systems has 
attracted much attention in the last two decades due to their applications in flight control [12], turbofan 
engines and wind turbine systems [13, 14]. For example, an application of the distributed control approach 
to control and synchronization of wind generation systems modeled as parameter varying systems has been 
presented in [14]. Even though a discrete time model was considered in [14], parameter varying system 
modeling was motivated by the fact that the wind energy source driving wind turbines exhibits time-varying 
nature. In order to integrate a wind generation system into a power grid in a grid-friendly manner, the total 
power output of the wind turbines must be regulated to conform to a constant output. Each turbine then 
serves as a node and is controlled with respect to its power output by turning the blade pitch angle, and the 
value of the blade pitch angle is propagated through the communication network using the leader-follower 
algorithm. Thanks to this and many other potential applications, the theory of cooperative control for 
parameter varying systems has been gaining attention in recent years [15, 16]. 

The main contribution of this paper concerns the leader-follower control of parameter varying intercon¬ 
nected systems. A related problem, from the synchronization viewpoint, has been studied in [18, 19]. In 
the first reference, the synchronization problem is solved for heterogeneous systems which depend on the 
parameters in an affine fashion, under the assumption that the network of LPV agents allows for an internal 
model for synchronization while the agents are decoupled from each other. That is, the agents interact over 
the control protocol only. In [19], while the linear dependency on the parameter is not required, the leader 
is assumed to be given and be completely decoupled from the agents. In contrast with these references, we 
assume the leader to be chosen from the group of parameter varying agents; it is interconnected to the rest 
of the network, and the linear dependency on the parameter is not required. 

The fact that the leader is interconnected with the followers makes it difficult to apply regular central¬ 
ized tracking techniques to the problem under consideration. Many tracking techniques assume that the 
reference trajectory is generated by an exosystem or is a priori known, and is independent of the followers. 
This is not the case in this paper. Alternatively, a centralized tracking controller can be obtained by solving 
a stabilization problem for a large-scale system comprised of subsystems describing dynamics of individual 
tracking errors. This can be done by applying one of the existing centralized gain-scheduling robust control 
techniques, e.g., using the technique from [21]. However, in general the centralized controller obtained this 
way will not have a desired information structure. Indeed, a controller obtained in such a way will not 
generally guarantee that the information from subsystems that are not observed by a node is not required 
for feedback. One way to enforce such an information constraint is to employ a block diagonal Lyapunov 
function. This is the solution approach undertaken in this paper and is its main difference from solutions 
which could be obtained in a centralized setup. 

Different from many papers that study synchronization of decoupled systems (including the above 
mentioned references [18, 19], also see [2, 3]), in the case of coupled systems, it is essential to distin¬ 
guish the network representing the existing interactions between the subsystems (including interactions 
with the leader) from the network which realizes communication and control. This leads us to consider 
a two-network structure in this paper. The rationale for this is twofold. Firstly, our aim is to construct 
synchronization protocols for all agents excluding the leader, as we wish to avoid perturbing the leader 
dynamics (other than through an unavoidable physical coupling). Hence, the communication graph of the 
network must be different from the interconnection graph. The second reason is that the interconnection 
graph describes dynamical couplings between subsystems and thus may be different from the communica¬ 
tion and control graph. 

Our main result is a sufficient condition for the design of a gain-scheduled leader-follower control 
protocol for parameter varying multi-agent systems with directed control network topology and linear 
uncertain couplings subject to norm-bounded constraints. The condition involves checking feasibility of 
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parameterized linear matrix inequalities (LMIs) at several operating points of the system. These LMIs 
serve as the basis for the design of a continuous (in the scheduling parameter) control protocol for coupled 
parameter varying systems, by interpolating consensus control protocols computed for those operating 
points; cf. [29, 19, 21]. Interpolation allows us to mitigate detrimental effects of transients arising when 
the system traverses from one operating condition to another. 

The remainder of the paper proceeds as follows. In Section 2, we formulate the leader follower control 
problem for parameter varying coupled multi-agent systems and give some preliminaries. The main results 
are given in Section 3. Section 4 gives an example which illustrates the theory presented in the paper. 
Finally, the conclusions are given in Section 5. 


2 Problem Formulation and Preliminaries 


2.1 Interconnection and communication graphs 

As stated in the introduction, in this paper we draw a distinction between the communication topology 
of the system used for control and the topology of interactions between subsystems. The example of the 
two-graph structure is shown in Fig. 1, where the edges of the interconnection graph are indicated by 
the solid lines and the edges of the communication graph are shown by the dashed lines. Both coupling 
and communication topologies are described in terms of directed graphs defined on the common node set 
V = {0,..., N}. Without loss of generality, node 0 will be assigned to be the leader of the network, 
while the nodes from the set Vq = N} will represent the followers. The coupling graph and the 

communication graph will be denoted as and C/'^, respectively. 

The edge sets of both graphs are subsets of the set V x V and consist of pairs of nodes. The pair (j, i) 
in each edge set denotes the directed edge which originates at node j and ends at node i. Edge (j, i) in the 
edge set of the directed coupling graph denoted £‘^, describes the fact that node i is influenced by node 
j through a directed interaction between nodes j and i. 

Also, f ° C V X V denotes an edge set of the communication graph Q‘^, consisting of ordered pairs of 
nodes. Each such edge indicates the information flow between the nodes, that is, {j, i) G if and only if 
node i obtains information from node j, which it can use for control. 

The adjacency matrix of the directed interconnection graph is denoted as its (i, j)-th entry is 
1 if and only if (j, i) G 8'^. The adjacency matrix of the directed communication graph is defined 
in the same manner. Since according to a standard convention we assume that both the coupling graph and 
communication graph have no self-loops, the diagonal entries of A'^ and A^^ are all equal to zero. 

To distinguish between the leader and the rest of the network, we define subgraphs Q^, Qq of the graphs 
0‘^, defined on the node set Vo = N}, with the edge sets = {(j, i) G £^ i, j G Vo} and 

£q = {(j, i) G £'^'. i,j G Vo}, respectively. Selecting the subgraphs Qq and Q§ induces the partition of the 
matrices A"^, 


A‘^ = 


'0 

d 

d 

'^0 J 


A^ 


'0 

0 ■ 

_g 

A^ 
^0 J 


where d = [di ... d^]', d = [di ... 9 = [gi ■ ■ ■ gN]' and Aq, Aq are the adjacency matrices of the 

subgraphs Qq and Qq, respectively. The zero row of the matrix A'^ reflects our assumption that the leader 
node does not receive the state information from other nodes of the network. However, the corresponding 
row of A^ may be nonzero since the leader node may be physically coupled with some of the followers. 

In this paper we are concerned with the case where only some of the followers receive the state infor¬ 
mation directly from the leader. We refer to such nodes i G Vo as pinned nodes; the corresponding entries 
of the adjacency matrix gi = 1, and = 0 if node i is not pinned. The matrix G = diagjpi} G 
is referred to as a pinning matrix. 

The Laplacian matrix of the subgraph Qq is defined as £§ = V — Aq, where V = diagjpi ,... ,pn} G 
j^NxN j-jjg in-degree matrix of Qq, i.e., the diagonal matrix, whose diagonal elements are the in-degrees 
of the corresponding nodes of the graph Qq, pi = X]_;=i foi' i = ■ ■ ■, N, where are the elements 

of the jth row of the matrix Aq . 

Einally, we give the definition and the notation for neighborhoods in the above graphs. 
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Figure 1: An example of the interconnection and communication graphs. 


Definition! Qq = is a communication graph, each node i € V represents a subsystem in the 

interconnected system and edge (j, i) € f ° means that node i obtain information from node j. 

Definition 2 is an interconnection graph, each node i & V represents a subsystem in the interconnected 
system and edge {j, i) G indicates that subsystem i is influenced by subsystem j through a directed 
interaction between subsystems i and j. 

Node j is called a neighbor of node i in the graph (or Q^, Q'^, Qq, respectively) if {j, i) G £^ 
or £q, respectively). The sets of neighbors of node i in the graphs and are denoted as 
Nf = {j\{j,i) G 1 ?'^}, and AT = {j\{j,i) G respectively. The sets of neighbors of node i in the 
subgraphs Gq and Qq are denoted as Sf = {j\{j, i) G £q} and S'f = {j\{j, i) G £q}, respectively. 

We conclude this subsection by specifying our standing requirements on the communication topology 
of the network which are assumed to hold throughout the paper. 

Assumption 1 The communication subgraph Gq contains a spanning tree, whose root node ir is the pinned 
node, i.e., gi^ > 0. 

Remark 1 Assumption 1 ensures that the information flows from the leader to its pinned followers, as well 
as between other followers [25, 28 ]. 


2.2 Notation 

Throughout the paper, the following notation will be used. 

5ft" and 5ft"^'" are a real Euclidean n-dimensional vector space and a space of real n x m matrices. 

r denotes an interval [pmin, Pmax] C 5ft. The symbols p, ps, ps, p , etc., will represent various points 
within the interval T, and T^, Fg will denote sets of points within F which will be defined later in the paper. 
Also, pit) is a scalar function, defined on [0, oo) and taking values in the interval F. This function will 
describe a time-varying scheduling parameter of the multi-agent system under consideration. 

Unless stated otherwise, the notations A(p), Yp, etc., will refer to matrices of appropriate dimension 
parameterized by p G F. 

For q G 5ft", diag{g} denotes the diagonal matrix with the entries of q as its diagonal elements. 

0 denotes the Kronnecker product of two matrices. 

Amax(’) and Aniin(') will denote the largest and the smallest eigenvalues of a real symmetric matrix. 

In is the N x N identity matrix, and Iat G 5ft^ is the vector whose all entries are equal to 1. When the 
dimension is clear from the context, the subscript N will be suppressed, and I will represent the identity 
matrix of an appropriate dimension. 

For two symmetric matrices X and Y of the same dimensions, X >Y, (X > Y) if and only if X — Y 
is positive semidefinite (positive definite). 

Consider the Laplacian matrix of the subgraph Gq, According to [28], under Assumption 1, the 
matrix £q -G G is nonsingular (also see [27]), and all the entries of the vector P = [t?i,... j'dAr]' = 
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{Cl + G)-H isf are positive. Also, the matrix 0 + G) + (£§ + G)'0 ^ is positive definite, where 

0 = diagjr^}. The following two constants will be used in the proofs of our results: 


O’ — (® ^(-^0 + G) + (£g + G)'Q , 

A = A„,ax ((/:§ +G)'0-2(£=+G)). (1) 


2.3 Problem Formulation 

The system under consideration consists of TV + 1 parameter-varying dynamical agents, coupled with 
their neighbors; the topology of interconnections is captured by the directed graph Q'^. Dynamics of the 
jth agent are described by the linear equation 

Xi = A{p{t))xi + BiUi + B 2 - Xi), i = 0,...,N, (2) 

j&Nf 


where Xi G 5ft" is the state of agent i, Ui G is the control input and p(t) is the time-varying parameter, 
which is available to all agents. It is assumed that p: [0, 00 ) —T C 5ft is a continuous function. Also in 
equation (2), Bi and B 2 are real matrices of appropriate dimensions, and A{p{t)) is the composition of a 
continuous matrix-valued function A(-): T — 5ft"^" and the function p{t) defined above. 

The functions ipij : [0, c») x 5ft" — >■ 5ft'", 

Lpij{t^X^ — ^ij(j‘)GijXj X G^ , 

describe how agent j influences the dynamics of agent i. Here G^ G 5ft'J“i^" and G 5ft'"^‘^‘^ are 
respectively constant and time-varying matrices. According to the above model, we focus on linear time- 
varying interactions between the subsystems whose strengths depend on the relative state of subsystem i 
with respect to subsystem j, Xj — Xi. Also, we will assume that the coefficients are uncertain, and 

satisfy the constraint 

< -f) Vf e [0,oo). (3) 

That is, we consider a class of norm-bounded uncertain interactions, which will be denoted by S. 

Since we have designated agent 0 to be the leader, and the rest of the agents to be controlled are to 
follow the leader, then according to (2) dynamics of the ith follower are described by the equation 

Xi =A{p{t))xi + BiUi + ^ <fij{t,Xj - Xi) + d^ipM{t,xo - Xi)^, i = ( 4 ) 


while the leader dynamics are given by the equation 

±0 = A{p{t))xo + B 2 ^ ipok{t,Xk - xo). (5) 

k: dk — 1 

Unlike agents i, i = 1,... N, the leader is not controlled, i.e., uq = 0. On the contrary, all other agents 
will be controlled to track node 0. 

Remark 2 In this paper, no a priori stability assumptions are made about the matrices A{p). This is in 
contrast to, for example, the synchronization problem considered in [22], where the dynamics of the leader 
were required to be stable or marginally stable. In the sequel, however, certain conditions, in an LMI 
form, will be imposed on the coefficients of the leader and followers dynamics. While we formally make 
no special provisions regarding these coefficients, except for the feasibility of these LMI conditions, for 
such LMIs to be feasible one can reasonably expect the overall system dynamics to have certain collective 
stabilizability and detectability properties; see Remark 3 below. 
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Define tracking performance associated with the control input w(*) = [ui(-)^ ■ • ■ ^Ar(*)T 


^ pOO 

J(m)=^/ {{xq - XiyQ{xo - Xi) + u'iRui) dt, (6) 

where Q = Q' > 0 and R = R' > 0 are given weighting matrices. In this paper we are concerned with 
the following problem. 

Problem 1 For each follower i,find gain schedule functions : F —>■ so that the control protocols 

of the form 

Ui = -Ki{p{t)) < ^ {xj - Xi) + gi{xo - Xi)> , ( 7 ) 

Uesf J 

ensure a bounded worst-case tracking performance, 

snpj{u) < const, (8) 

where const denotes a constant which can depend on a;o(0) and Xi(ff). 


3 The Main Results 

In this section, we first revisit the leader follower tracking control problem for multi-agent systems 
with fixed parameters considered in [ 17] to obtain a sufficient condition for the design of a control protocol 
for a more general class of systems involving coupling between the leader and the agents. This sufficient 
condition will then be applied to develop an interpolation technique to obtain a continuous gain scheduling 
tracking control protocol for the system (2). The interpolation technique is the main result of this paper. 

3.1 Leader follower control for fixed parameter systems 

We now revisit the results of [17]. This revision is prompted by a more general structure of the system 
(2) which allows for physical connections between the leader and the followers. As a result, the form of 
the control protocol is somewhat different here in that the resulting controller gains depend on i. 

Consider a fixed-parameter version of the system (2) described by the equation 

Xi =A{p)xi-^ j3ii{t,Xi)BiUi-^ B 2 Pij{t,Xj - Xi), i = 0,...,N, (9) 

jGNf 


where p S F is fixed, and /3 is a positive constant. Compared to (2), the system (9) includes an additional 
uncertainty element Xi), which satisfies the following constraint 

Uo(,t,x) - <\\x - yW^, Va;,yG5R”. (10) 

In the sequel, we will show that when the difference A{p{t)) — A{p) is sufficiently small, the system (2) 
can be represented as the system (9) subject to (10). 

We now derive a distributed protocol of the form (7) under which the fixed-parameter uncertain system 
(9) satisfies the performance requirement (8). 

Define the leader tracking error vectors as = xq — Xi, i = 1,..., iV. Dynamics of the variable 
satisfy the equation 


Ci =A{p)ei - Biu^ I3ii{t) - B 2 ^ ipokit, Ck) 

k : !ij5 = l 

- B 2 ^ {(f^j{t,ei) - - B2di(f^o{t,ei), 


( 11 ) 
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where ^i(t) = ^o{t, xo{t)) — Xi{t)). It follows from (10) that |||i|p < ||ei|p for all f > 0. 

For node i of the subgraph Qq, introduce matrices Ci = Ci = J', where 

ji) • ■ •) jKi are the elements of the neighborhood set Sf, and ri,... ,r^^ are the nodes with the property 
(j, rj S £q', Ki and Xi are, respectively, the in-degree and the out-degree of node i in the graph Qq. Also, 
let Q = {a'^ /X)Q, where a, A are the constants defined in (1). 

In order to formulate the extension of Theorem 1 in [17], with each node i,i = 1,..., N,we associate 
a collection of positive constants , Hij, j G Sf, tti, i^io (only for those nodes i for which di = 1), and 
fiQi (only for those nodes i for which di = 1). Also, let F = F' > 0 be an n x n matrix. Using these 
constants and the matrix, for each node i introduce a matrix Hi defined depending on di, di as follows. 

Case 1. di ^ 0 and di ^ 0. For each such node i, define the matrix Hi as 


where 




YQ^!'^ 

YC'i 

YC'i 

Y 

yc',a 

o 

o 


Qdl'^Y 

-I 

0 

0 

0 

0 

0 


CiY 

0 


0 

0 

0 

0 

n* = 

CiY 

0 

0 

-Hi 

0 

0 

0 


Y 

0 

0 

0 

-Ai 

0 

0 


CioY 

0 

0 

0 

0 

-^7 

0 


CoiY 

0 

0 

0 

0 

0 

-jy^l 


=diag[—/, j G Sf], n, = diag[—/, j : i G Sf], 

^ij 

Z, =A{p)Y + YA{p)' - BiR-^B[ + -/J^J 


(—I—) -I—I- X/ —)^2.b2- 




fc: dk = l 


A. 



( 12 ) 


(13) 


Case 2. di = 0 and di ^ 0. For each such node i, the constant ViQ is not defined. Accordingly, we 
define the matrix 11^ by removing the second last column and row from the matrix in (12): 


Z^ 

FQI/2 

FC' 

FC' 

F 

YC',, 1 

Q^/'^Y 

-I 

0 

0 

0 

0 

CiY 

0 


0 

0 

0 

CiY 

0 

0 


0 

0 

Y 

0 

0 

0 

A. 

0 

Co^F 

0 

0 

0 

0 

WMoi ^ -1 


where the matrix Zi is modified to be 

Z,=A{p)Y + YA{p)'-B^R-^B[ + -l3‘^I+{Y.i— + —)+ Y. —)B2Bf (15) 

TT, jest = 


Case 3. di 0 and di = 0. For each such node the constant poi is not defined, hence the corresponding 
matrix 11^ will be defined by removing the last column and row from the matrix in (12): 



r 

YQ^C 

FC' 

FC' 

F 

1- 

o 

>-< 


qACy 

-I 

0 

0 

0 

0 

n. = 

CiY 

0 


0 

0 

0 

CiY 

0 

0 

-U, 

0 

0 


Y 

0 

0 

0 

-A, 

0 


_ C,oY 

0 

0 

0 

0 

—1-/ 

J 


The matrix Zi is the same as in (13). 
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Case 4. di = 0 and di = 0. In this case, both i^io and /ioi are not dehned, and the corresponding 
matrix 11^ is defined by removing two last columns and rows from the matrix in (12): 


n. = 


z, rOi/2 yc' YC' Y ' 

Qdf^Y -/ 0 0 0 

CjY 0 0 0 

C^Y 0 0 -Oi 0 

Y 0 0 0 -A, _ 


The matrix Zi is the same as in (15). 

The following theorem is an extension of Theorem 1 in [17]. 


(17) 


Theorem 1 Under Assumption 1, let a matrix Y = Y' > 0, Y G constants Vij > 0, pij > 0, 

TTi > 0, j S Sf, i = ... ,N, and constants ViQ > 0 (for those i with di f 0), poi > 0 (for those i with 

di 0) exist such that the following LMIs are satisfied simultaneously 


Hi < 0, * = 1,..., A^. 


(18) 


Then the control protocol (7) with 


Kfip) = -{d,a)-^R-^B[Y-^ (19) 

solves the leader follower tracking control problem for the fixed parameter system (9). Furthermore, this 
protocol guarantees the following performance bound 

A ^ 

sup J(m) < — ^(a:o(0) - a;j(0))'y ^(a:o(0) - 2 :^( 0 )) (20) 

“ i=l 

for all uncertainties ^ifor which (10) holds. 

Remarks The feasibility of the LMIs (18) in Theorem 1 is a sufficient condition for the existence of a 
control scheme to guarantee the performance bound (20). The LMIs are numerically tractable and can be 
solved using the existing software. Additionally, using the standard tools of the Hoo control theory, such as 
the Strict Bounded Real Lemma and the associated Riccati equation, the feasibility of the LMIs (18) can be 
related to collective stabilizability and/or detectability properties of the coefficients of the interconnected 
system (9). Since this system consists of identical agents, one can expect that each agent must necessar¬ 
ily have corresponding stabilizability and detectability properties for these collective properties to hold; 
cf [20 f 

Proof: Using the Schur complement, each LMI (18) can be transformed into the following Riccati 
inequality: 

A(p)y +rA(p)'-Si+ ( ^ (— + —)+ ^ 

+ y((3+ ^ ^ VijC[jCij-\- 'y ^ p,jiCOCjiTTii^Y 

+ ( — B2B2 + Vi^YCl^CufY 
\Vio 

(this term is present only if di = 1) 

+ NpoiYC'ofioiY < 0 (21) 

(this term is present only if di = V). 

Note that the last and the second last lines in the Riccati inequality (21) are present only for those nodes i 
for which di f 0 and/or di f 0, respectively. 
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After pre- and post-multiplying (21) by y and then using the expression (19) for Ki{p) in the 
resulting inequality, we obtain 

^ + {A{p) + a'diBiKi{p))'Y VijC[jCij + pjiC'^^Cji 

jesf j-.iesf 

+ y-^f( y (— + —)+ y —)B2B!, + -P^i)y-^+Y-^BiR-^B[Y-^+Q + ttJ 

+ ( - Y ^ B 2 B 2 Y ^ + I'ioCioCiO^ 

/ 

(this term is present only if di = 1) 

+ Npo^CUCo^ < 0 (22) 

(this term is present only if di = 1). 

Define e = [e'l,..., e^]' and consider the following Lyapunov function candidate for the intercon¬ 
nected system consisting of the subsystems (11): 


N 


y(e) = ye'y- 


( 23 ) 


Then 


dV{e) 

dt 


N 

= y 2e'y"^(A(p)e, -f BiA:,(p)(y (ci - ej) + giCi)^ 


2=1 


j^Sf 


N 


N 


N 


- 2 y y e'r ^B2<Pofc(f,efc)-t-2y e'y Ci) - 2 y djC-y ^B2(pio{t,ei) 

N N 

- 2y y e'y-ii? 2 <^*,(t,e,) + 2y y e'y-iB 2 ^.,(f,e,). (24) 

i=i jest i=i jest 


Note the following inequality: 


N 

y 2 e-y^BiA:,(p)( y (Ci - Cj) + Pie*) 

i=l j^S'r 

= -2e\Q{Cl + G)® {Y-^Bi{aR)-^B[Y-^))e 

= -y'(( 0 (£S + G) + {Cl + G)' 0 ) ® 4)2/ 

< -2ay'{lN O 4 )?/ 

= -2ae'{iN 0 y^Bi(cri?)-iB;y-i)e 

N 

= -2 y e'y i^iy 4 i, ( 25 ) 

2=1 


where p = (Jat (8) (cri?) i)e. 
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From (24) and (25), one has 


dV{e) 

dt 


N 


N 


i —1 i —1 

N N 

-2EE e'y ^B2ipij{t,ei) - 2E E e-F ^B 2 (pQk{t,ek) 

i=l jeSf i=lfe:dfc = l 

N N 

- 2'^d^e'^V~^B2<p^o{t,ei) + 2^^ ^ e[Y~^B2(pij{t,ej). 
i=i i=i jesf 


Using the Riccati inequality (22), it follows from (26) that 


(26) 


dU(e) 

dt 


N 


j(^sf 


idSt 


+ ^,/ + F-i((^(—+ —)+ ^ J-)B2B'^ + -fi^l)Y-^ 


N 


N 


+ E! ^B 2 B 2 Y ^ Y ViQC'^QCio^ei + N E! e^/roiCgjCoiei 

di = l i: di = l 

N N N 

Y 2 Y,e',Y-^PUt,ei )-2 EE e-F ^B2(pij{t,ei) - 2'^die'iY ^B2ipio{t,et) 

i^l i^l j£Sf 


iV 


iV 


-2E E + 2E E ^B 2 (pij{t,ej). 

i=lfc:dfc=l i=lieSf 


(27) 


Using the following identities. 


N N 

E E E d^e'C'^C.^e,, 

i=ljgSf *=li:ieS|’ 

N N 

^ E ^i/^oiC'oiC'oie- =A^ E efcA^ofcC'g^C'o/cefe 

2: (Ji = l fe : (ife = l 

N 

'^die[Y~^B2Pio{t,ei) = E e'F“^S2(/7io(i, e^), 


2 = 1 


2: di = l 
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and completing the squares, one has 


dV{e) 

dt 


< - 


N 


1 








N 


- EII ®*)ii' + E e*)ii' - iie^ii') 


N 


z=l 

N 


i=l 


N 


-EE \\—=B2Y ^ei + y/u~ipij{t,ei )ii“ + EE>'«(ii*’«(‘.«.)ii"-iic«e,ii“) 

i^l J^sf V i^i j^sf 

N N 

-EE II ! — —^2^ 's/ ) II EE P'ij)(||</5y (t,ej)|| IlC'ijejIl ) 

i=l jgs|’ 


i=l jGSt 


1 


di = l 


N 

-E E II . - ^2^ + •v//tofeV^ofe(i, efc)||^ + E A^ofc(||v^ofc(i5 efe)||^ — IjCofeefcll^)- 

fclfc :4 = l fe :4 = l 

(28) 


^ei + ei)||^ + ^ ^'io(||v5zo(^, ei)|p - HCioeilH 

1 


According to the norm-bounded condition (3), from (28) we have 








— 1 D' 1 




q)< 


Since y(e(f)) > 0, then (29) implies 

+Q)eidf < y(e(0)). 


dt. 


(29) 


JV 

S/ 

z=l -^0 


(30) 


The expression on the right hand side of the above inequality is independent of t. Letting t ^ oo leads 


to 


e',(Y-^BiR-^B[Y-^+Qy,dt<V{e{0)). (31) 

Using (6) and (7), we have 

^e-Qei + u'^Rui^dt 

(e'(/jv 0 Q)e + e'{m + + G) (g) ^Y-^B^R-^B[Y-^)eyt 

/•oo X \ 

; J y{lN(gQ)e + e{lN(g-^Y-^BiR-^B[Y-^)eyt 


N 

Jiu)=Y. / 

i=l -^0 


=v r^j^Y- 


i=Eo 


e',(A^Y-^BiR-^B[Y-^ + g)e. 


dt. 


(32) 


Since Q = -^Q, it follows from (30) that 

Af / 


\ N poo \ ^ 

J{n)<—Y. / e'(y-ii7ii?-ii3;y-i + Q)e, < ^^e'(0)y-ie,(0). 

^ i=l •'0 ^ i=l 


(33) 
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It implies that the control protocol (7) with Ki{p) defined in (19) solves leader following tracking control 
problem, and also guarantees the performance bound (20). □ 

Theorem 1 can be applied to obtain a leader-follower tracking protocol (7) for the system (2) if param¬ 
eter variations of systems are sufficiently small. 

Suppose there exists po G F and /3 > 0 such that 

{A{p{t)) - A{po))'{A{p{t)) - A{po)) < P'^I, V f G [0, oo). (34) 

and define Xi) := ^[A{p{t)) — A{pQ)\xi. This allows us to regard small variations of the matrix A(-) 
as perturbations. The fixed-parameter system (9) with p = po captures this type of perturbations. Then the 
following result follows from Theorem 1 . 

Corollary 1 Under Assumption 1, if the LMIs (18) with p = pQ and Y = Ypg are satisfied simultaneously, 
then the control protocol (7) with 

Kfipo) = -ii},a)-^R-^B[Y-^ (35) 

solves the leader following tracking control problem for the parameter varying system (2) under small 
variations of pit) for which condition (34) holds for all f > 0. Furthermore, this protocol guarantees the 
following performance bound 

A ^ 

supj(u)< — ^(a:o(0)-a;i(0))'y'^^^(xo(0)-a:i(0)). (36) 

3.2 Design of a continuous protocol schedule 

The result of Corollary 1 only holds under assumption that variations of the matrix A(p( )) are suf¬ 
ficiently small to satisfy (34). In general, it may be difficult to satisfy (34) using a single p^, or the 
cotTesponding LMIs of Corollary 1 may not be feasible. In order to address this situation, we propose a 
gain scheduling approach. 

Consider a set of design points F^ := {ps,s = C F and a collection of positive constants 

/3s chosen so that for any p G F there exists at least one point ps with the property 

(A(p) - Al(ps))'(A(p) - A(ps)) < fill (37) 

and that the LMIs (18) with p = Ps are feasible. 

Let Us be the largest connected neighborhood of the design point ps G F^ such that (37) holds if 
p = p{t) G Us- This allows a protocol (7) to be scheduled for each value p on the trajectory p(f), by 
associating with every p{t) the protocol computed using Theorem 1 for one of the indexes s G {s: p(t) G 
Us}. However, when applied to the parameter-varying system (2), the gains of such a protocol may become 
discontinuous at the time instant when p{t) is switching between different sets Us- To overcome this issue, 
the continuous interpolation technique proposed in [29, 21] is used in this paper to obtain a continuous 
consensus control protocol; also see [19]. 

Consider an arbitrary fixed p G F, and the collection of constants fis and grid points F^. Select s such 
that p G Us, and let , z/y^p^, p^-^p^, Yp^ , {uio^p^ , poi,p.), j G Sf , i = 1, ..., N, be a feasible solution 
to the LMIs (18). Recall that Vio,p„ and poi,p, are only defined for those nodes i with di f 0 and di f 0, 
respectively. It is straightforward to show that this collection of positive constants and the matrix is also a 
feasible solution to the following reduced coupled LMIs 

T, <0, r = l,...,/V, (38) 

where the matrix is defined as follows. For those nodes i, for which di Q and ^ 0, is defined 
by removing the fifth column and row from the matrix H^ defined in (12) and replacing the matrix Zi in 
(12) with the following matrix: 

Zi = A{p)Y + YA[p)' — BiR ^B'i+ ( (-1-) H-h \B 2 B 2 . 

Stithy M.. ^^o 
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where the matrices <i>i and fli are as defined previously. That is, the new matrix Zi is obtained from the 
matrix in (13) by subtracting the term This results in the matrix defined as 



YQ^/"^ YC' YC' 

YC',, 

YC(,, 1 

QY'^Y 

-I 

0 

0 

0 

0 

aY 

0 


0 

0 

0 

C,Y 

0 

0 


0 

0 

C^oY 

0 

0 

0 

— 

0 

Co^Y 

0 

0 

0 

0 

Npoi ^ \ 


For three other cases, the matrix is defined in the same fashion. First, the fifth column and row are 
removed from the matrix 11^ defined in (14), (16), or (17), respectively. Next, the matrix Zi is redefined by 
subtracting from the corresponding matrix in (13) or (15), as appropriate. 

Now consider the uncertain fixed parameter system (11), and assume p G UsnUs+i- Then we conclude 




that both collections , fJ'ij^p ^, Yp ^, Ir'io.ps; Fot.ps) ± p^.,.!, 

jG5r,i = i,...,iv, are feasible solutions to the coupled LMIs (38). 

This allows us to construct interpolated feasible solutions to (38) as follows. 

For a 7 G [0,1], define 


t'iO, 


Ps + l ’ 


fJ'0i,Ps+i )’ 


>"7 = 7yp.+(l-7)i"p.+i, (40) 

“ 7)i^*7p,+J“\ (41) 

= [7M,7p, + (1 - 7)li[;]-,p,+J“^- (42) 

Also, for nodes i with di = 1, define 

= bi^rolp, + (1 “ l)^io]ps+i\~^ ^ ( 43 ) 

and for nodes i with di = 1, define 

= [7/io*!p. + (1 - (44) 


Lcitimtl 1 Given i f-^ij^ps jYp^ ^ |^t0,Ps5M0z,ps) ^^d f^ij^p^^nYp^^^^ \^iO,pE+i: i^ 0 i,pE+i')> j ^ 

Sf, satisfying the LMI (38). Then {I'ij^-y, Yj, |i^io, 7 i Mot, 7 ) olso satisfies the LMI (38). 

Proof: The statement of this lemma follows from the observation that the inequality (38) is linear with 
respect to the variables Y and (—-—, —-— and —-—, —-—) where these latter variables appear. □ 
Using the above lemma, we now define a collection of interpolated gains for the control protocol (7), 
as follows. Suppose the collection of positive constants /3s, fis+i and the grid points Ti have the following 
properties: If p G [ps,Ps+i], then 

(A{p) - A(p,))'(A{p) - A(ps)) < fill, Ps<p<Ps, (45) 

(A(p) - A(ps+i))'(A(p) - A(ps+i)) < fil+il, p^_^^ < p < Ps+i, (46) 
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where ps < p^^^ < Ps < Ps+i- Define 7 = '){p) = and let 


s + 1 


'Yps, pe [ps,p^^J, 
Yp = {y^, pe [p^^^,Ps], 

^s + ll P^ (PsjPs+i]i 


'^3,P 


Also, for nodes i such that di = 1, define 


''’■J,Ps ) 


^ii,Ps+ii PG (PsiPs+l]j 


Pij,Pa1 p£ [^**’£ 8 + 1 )’ 

“ S Pij,!-: P ^ [^s+i’^s]’ 

^Pij,Ps+n P^ (Ps;Ps+l]i 


(47) 


(48) 


(49) 


^iO,Ps^ PS [PsiP 


■s+l^ 

l^iO.p = ^ ^'i0,7j P G 1^8+1’ 

,^i0,ps+ii PS (PsiP8+l]- 

Likewise, for nodes i such that di = 1, define 


(50) 


P0i,PsJ PS [/^^’^8+l^’ 
p0i,p = ^ p0i,7) P G I^s+l’^®]’ 

,p0i,Ps+i) PG (Ps)P8+l]- 


Next, define the gain for the control protocol (7) 


7fz(p) = 


r-l 


(51) 


(52) 


The function Ki is a continuous function on T, since Yp > 0 for all 7 G [0, 1]. Let Lq be a set 
consisting of all the corner points p^^j^,P8 G T. Without loss of generality, it is assumed that the set 
{f > 0: p{t) G Lq} has zero Lebesgue measure. 

The following theorem is the main result of this paper. 

Theorem 2 Under Assumption 1, suppose that the time-varying parameter p(-) of the uncertain linear 
system (2) satisfies the condition 


j-j ^ 

sup|p(t)| < -g, Q= sup 

t Q pGr\ro 


dy-i 
I—^1 

' dp ' 


(53) 


where p = min Amin (l^p ^BiR ^B'^Yp ^ + Q) and q G [0,1) is a constant. 

Then the control protocol (7) with the gain schedule Kif) of the form (52) solves the leader follow¬ 
ing tracking control Problem 1 for the system (2). Furthermore, this protocol guarantees the following 
performance bound 


A ^ 

supJ^(m) < ——^^(a;o(0) -a:i(0))'(rp(o))"^(a;o(0) -X, 


( 0 )). 


(54) 


Proof: Since the matrix Yp ^ is continuous and piecewise differentiable except at p G Lq, it follows 
from [29] that, given any e > 0, there exists a continuous differentiable matrix function Xp defined on T, 
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and a constant <: > 0 for any corner points Pc S Fq such that 


sup||Xp-r II < £, pe + 


per 


^p = ^p^ pi{Pc-‘;,Pc + ^), 


I II II ^^p 


-1 


= max 


sup ||— 7^11 < sup ,1 , ,1 —„ , 

per dp per\ro dp dp 


P I 

sup II- ' 


= _ llr^ ^{Yp^ - Yp^+JF^ 


-ll 


s=l,...,M 7e[0,l] 


(55) 

(56) 


(57) 


Note that the approximating matrix Xp can be chosen symmetric, since Y~^ is symmetric. Also by select¬ 
ing a sufficiently small £ > 0, a positive definite matrix Xp can be selected for all p S F. 

Consider the closed loop large-scale interconnected system describing tracking error dynamics of the 
system (4) and (5) 


ti = A{p{t))ei + BiK^{p){'^ {et - ej) + giet) 

- 82 '^ tpij{t,ei-Cj) - B 2 ^ (pok{t,ek) - B 2 d^pio{t,ei). (58) 

jeSf k: dfc = l 

Let the following Lyapunov function candidate for this system be chosen 

N 

l/(e) = ^e'Xpe,. (59) 

i=l 


We have 


dV{e) 

dt 


N 

■-'^2e^Xp(^A{p)ei + BiKi{p){ ^ (ci - ej) + giei)J 
t=i jeS^ 

N N 

-2E E e^XpB2(pok{t, eu) -f 2EE e[XpB2ipij{t,ej) 

*=1 k: dk = l »=1 ieSf 

N N N 

2^^ ^ ^ C^XpB2^ij {t ^ €i') 2 ^ ^ djCj^XpB2Pm{t^ Cj) ^ ^ GjXpCj . 

i=l j^sf i=l i=l 


(60) 


N 


Next, a bound on 2e'Xpi?iiFi(p)( (^i 


jesf 


Sion 


Cj) + giCi) is obtained. First we transform this expres- 


N 

2e[XpBiKi{p){Y (et - Sj) -f g^e^) 

i=l jGS‘ 

= -2e'(0(£= + G) ® iXpB,{aR)-^B[YY))e 
= -2e'{em + G) 0 {XpBi{aR)-^B[Xp))e 

+ 2e'(em + G)(^{XpBi{aR)-^B[{Xp-Yp-^))y. (61) 
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Let J = {aR) and consider 


+ G)'0) 0 JXp - 0 JiXp - r-i))' 

X ((^/i(£S + eye) ® jXp - -^In ® J{Xp - r-i)) 

= £(0(£§ + G)(£g + G)'0) ® XpJ'JXp + ^In ® {Xp - Yp-^)J'J{Xp - Yp-^) 

- (0(/:s + G)) ® XpJ'J{Xp - r-1) - ((£g + G)'0) 0 {Xp - Y-^)J'JXp > 0. 

It follows from (61) and (62) that 

N 

^ 2e[XpBiK,{p){ ^ (e, - e,) + gtei) 

i=i jeSf 

< -2e'(0(£g + G) 0 {XpJ'JXp))e 

+ e'(e{Qm + G)(£S + G)'0) 0 XpJ'JXp + 0 {Xp - Y-^)J'J{Xp - y-i))e 

= -y'{H 0 Ip)y + e'(^e{e{Q + G){C^o + G)'0) 0 XpJ'JXp 
+ 1/jv 0 {Xp - Y-^)J'J{Xp - y-i))e 

N N 

< -2ay'{lN 0 /p)y + ^eCll IlGf + Y.^\\Jr\\ey\^ 

N N 

= -2<je'{lN 0 XpB^{aR)-^B'^Xp)e + ^ eCl^^pf l|e,f + ^ e||^f ||e.|P 

i—1 i—1 

N N N 

= -2Y,e{XpB,R-^B',Xpe. + ^ eCl^^pf l|e.f + ^ 

i—1 i—1 i—1 

where y = {In 0 {aR)~^GB'^Xp)e and ( = Amin(0(£S + G){Cq + G)'0). 

Substituting (63) into (60), we have 


dV{e) 

dt 


N 


N 


< 


J2HXp{A{p)-B,R-^B',Xp)e, + 2Y, E e',XpB2^^j{t,e,) 

2-1 Z-1 j^Sf 

N N 

-2E E (f; 2 ^ ^ diCj^XpB2^iO (t^ 62) 

k: 4-1 

AT AT 

- 2 EE ^i^p^2^ij id-i ^i) Ee'A.e, 

2^12-1 


N AT 

+ E£Cl|J^pfl|e,f+ E^II'^II'II®*II'- 

2-1 2-1 


(62) 


(63) 


(64) 
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Consider the expression ^ 2e[Xp(^A{p) — BiR ^B[X^ei 

N 

- B^R-^B[Xpy, 

i^l 

N 

^ ^ 2e'Xp(^A{p) - B,R-^B[{Xp - Y-^ + y-^)) 

i^l 

N N 

< J2 2e:^P [a{p) - B,R-^B[Y-^) e, + ^ eae[ 

i^l i^l 

N 


-.^XpfJXpC^ 


+ ^ -ae[{Xp - Y-^)J'J{Xp - Y-^)e, 
i^l ^ 

N N N 

< ^2e-Xp^A(p) + adiBiKi{p)^ei + ^ecrU JXp||^||ej||^ + ^ecrU J||^||ei|| 
i—1 i—1 i—1 

Substituting (65) into (64), we obtain 


(65) 


<y 2e'Xp(^A(p)+cn9,Siit:,(p))e* + 2^ ^ e[XpB2(ptjit, Cj) 

i=l i=ljesf 

N N N 

-2E E ^iXpB2^Qk(i-) ^ ^ ^ di€j^XpB2^io(t^ Cj') -h ^ ^ ^jXpCj 

■i — 1 7^ . J n -i — 1 0 — 1 


( 66 ) 


^—1 k: dk — 1 
N 


N 


2 = 1 
N 


iV iV iV 

- 2^ ^ e'Xpi? 2 ^,,(t, e.) + E£(C + a)|| JXpf I|e,f + E ^(1 + ^)llH^.f • 
i=l j6Sf *=1 j=l 

Next, we turn our attention to the LMIs (38). Using the Schur complement, each LMI (38) can be 
transformed into the following Riccati inequality 

A{p)Yp + YpA{py-B,R-^B[ + {'£{— + ^)+ ^ -2—)B2B'^ 


+ YpiQ + E/ + E/ ^i3,P^ij^i3)^P 

3-^&Sf jGSf 


+ ( —^ 2^2 + ^^o,pYpCioCMYp 

\ ^iO,p 

(this term is present only if di = 1) 

+ Npo^,pYpC'o,Co^Yp < 0 

(this term is present only if di = 1). 

Note that the last two terms only appear in the Riccati inequality (67) for those nodes i for which = 1 


(67) 
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and/or di = 1. After pre- and post-multiplying (67) by Yp ^ and substituting (52) into it, we obtain 


Y-\A{p) + d,aBiK,{p)) + {A{p) + d,aBiK,{p))'Y-^ + Yp^ B^R-^ B[Y-^ + Q 

+ ( E + 77^) + E 77^) + E + E 

+ (—Y-^B2B'^Yp^ + v^o.pCioC^o'] 

\ ^iO,p / 

(this term is present only if di = 1) 

+ -^/tOi,pC'oi^Oi < 0 

(this term is present only if di = 1). 


(68) 


Since the set F is compact and coefficients of the Riccati inequality (68) are continuous in p, then 
provided e in (55) is sufficiently small, replacing Y~^ with Xp in (68), except the term Y~^ BiR~^ B[Y~^, 
preserves the strict inequality: 

Xp{A{p) + d,aBiK,{p)) + {A{p) + d,aBiK,{p))'Xp + Yp^ B^R-^ B[Yp^ + Q 

+ ( E ( 77 ^ + 77^) + E ^)XpB2B',Xp + E + E 


/esf 


jesr fe: 4=1 

+ ( - XpB2B2Xp + i'io,pC'iQCio\ 

\ ^iO,p / 

(this term is present only if di = 1) 

+ Xpoi^pC'^iCoi < 0 

(this term is present only if di = 1). 

Using the Riccati inequality (69), and the following identities 

N N 

j=ijesf *=ii:ieSf 


3- i&Sf 


(69) 


N 


N 


X ^ ^ 6j/roi,pC'oi6^0i6i =-^ ^ ^ Sfc/^Ofc.pC'ofc^^Ofcefc, 

2: = l fe : = l 

N 

'^die'iXpB2(pio[t,ei) = E e'^XpB2ipio{t,ei), 


i—1 i: di — 1 

and completing the squares, it follows from (66) that 

N „t 


pt JT// \ pi ^ pi 

/ e'(F-ii?ii?-iS7y-i+Q)e.df + E / e'(E) 

N N 

+ E E '^^3,pi\\‘Pi3it>ei)f - WCijCiW^) + Y^ E - IlCy-e^-f) 

i=i/esf i=i jesf 

+ E ^»o,p(||'^io(/,e*)|P - IlCioe^lP)+/V ^ /iofc,p(||/?ofc(/, efe)|p - IlCopefeU^) 

i: cii = l fe: 4 = 1 

N i-t 

+ E/ e((C + a)l|J^pf + (l + ^)l|Jf)l|eifd/. (70) 

7o 


idt 
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Furthermore, using the norm-bounded condition (3), we obtain 


-^2j +' 


N 


cT)\\jr)\\e,rdt + 


e'{x,y 


It follows from (53) that Y^ ^ satisfies the condition 


||17i < qXminiYp-^BiR-^B[Y-^ + Q). 
Together with (57), this inequality yields 


^dt. (71) 


(72) 


Since 


we have 


llTfpll < qX^iniY-^B^R-^B[Y-^+Q). 


X^i^iY-^BiR-^B[Y-^ + Q)I < Y-^BiR-^B[Y-^ + Q, 


V{e{t)) - V{em < - E / e' (y-^B,R-^B[Y-^ + q) < 


idt 


i=i-^o 

N 


+ 9 E/ e',{Y-^BiR-^B[Y-^+Q)edt 


2=1 

+ E / ^((C + 0')ll'^2^plP + (1 + ' 


cj)\\J\\^)\\ei\\^dt. 


Since l^(e(f)) > 0, (74) implies 

N 


(1 - g) E / e’,{Y-^B,R-^B[Y-^ + Q)e,dt 

i=i -^0 
N i-t 

-E I ^((C + a)||JXpf+ (l + a)|lJf)||e.fdf<y(e(0)). 


We now choose £ > 0 to be sufficiently small to ensure that 


(73) 


(74) 


(75) 


£1 = (l-g)77-£ |^(C + cr)(max||rp ||+£) + (1 + cr) 

is positive. Such an £ > 0 exists since at £ = 0, £i = (1 — q)rj > 0, and as function of £, £i is continuous 
at £ = 0. Then with this £, 

0 < £i < (1 - g)77 - £ ((C + a)\\JXpf + (1 + cr)|| Jf) , 
and it follows from (75) that 

N pt . N 


E/ \\eifdt<^'^{\\{Yp^o)) i+£)||e, 
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The above inequality holds for alH > 0 and the right-hand side is independent of t, therefore we conclude 
that lim(_>.oo /q ||ei||^df exists and is finite. This allows us to let f —>■ oo in (75) to obtain 


poo 

(1 - g) ^ + Q)edt 

N i-oc N 

<E / e((C + a)||JX,f+ (l + a)||jf)||e,fdf + ^e.(0)%(o)e.(0). (76) 

Jo 


Note that the left hand side of (76) is independent of e. Then we can let £ —?> 0 in (76). Since 
^p{o) (^p(o))~^ as £ —>• 0, this leads to 

N .oo _ 1 ^ 

Y, e[{YYB^R-^B[YY+Q)e^dt<- -^ eK0)(r,(o))-ie,(0). (77) 

i=i -^0 ^ ® i=l 


Using (6) and (7), we have 


^ POO C 

J{u) < V / e^d—YYBiR-^B[YY + Q)e^dt 
. 1 7o ^ 


2=1 ‘ 


e',{YYBiR-^B[YY + Q)edt 


< 


X 


N 


s2 E(^o(0) - a;*(0))'(yp(o)) ^(a:o(0) - a;i(0)). 

^ i—1 


(78) 


It implies that the control protocol (7) with Ki{-) of the form (52) solves leader following tracking 
control Problem 1 for the system (2), and also guarantees the performance bound (54). □ 


Remark 4 It should be noted that the proposed solution depends on the global information on the system 
topology. Specifically, the constants 9i and a are determined by the communication topology (but do not 
depend on the interconnection topology). On the other hand, the LMIs (18) are setup using the knowledge 
of interconnection topology. 


4 Example 

To illustrate the proposed method, consider a mass-spring-damper system in Fig. 2. The system consists 
of 21 identical masses which are coupled by different springs and dampers. Each mass is also connected 
to the wall with a spring. The dynamics of the coupled system are governed by the following equations 


mao —ki^i{ai — ao) + fci,2(<5i ~ <5o) + ^ 0 , 1 ( 0^20 ~ cto) + fco,2(<520 ~ Q^o) ~ kao, 
mai — ki^\(ai—\ a^) -t- k.i 2 (,d^i—i 

+ h+i^i(ai+i - a^) + ki+i^2iat+i - at) - kat + Ut, z = l,...,19, (79) 

ma2o =^ 20 , 1(^19 — q;2o) + feo,2(di9 — Q!2 o) + ko,i{ao — a2o) + ^ 0 , 2(09 — d2o) — ^020 + U 20 , 

where ai, ai are the displacement and the velocity of mass i, k and 1, i = 0,..., 20 are the spring 
constants, ki^2, J = 0,..., 20 are the damper coefficients, and m is the mass of each block. 

Choosing the state vectors as Xi = (ai, ai), i = 0,..., 20, each equation in (79) can be written in the 
form of (4), (5), where 


r 0 n 

, i?i = 

■ 0 ■ 

, B2 = 

■ 0 ■ 

1 

0 


J_ 


. m . 


. m . 


(fii^i-i (z) = (pi-i^i{z) = [kip kip]z, (z = 1,..., 20) 
Vo, 2 o{z) = V> 2 Q,o(z) = [ko,i ko, 2 ]z, z e R^. 
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Figure 2: Mass-spring-damper system. 



Figure 3; Undirected coupling graph. 


To verify the norm-bounded condition, let us assume that 

[ki,l ki^2\ = [ko.l fco,2] = Ao,2oCo,20j 

where Ci^i-i = [ki^i ki^ 2 \, Co ,20 = [^o.i ^ 0 , 2 ], are given matrices, and Ao ,20 S [0,1] are scalars. 

Then all the couplings between the agents satisfy the norm bound condition (3). In this example, for 
simplicity we let 1 = 2 = k for alH = 0,..., 20. 

The structure of the system (79) suggests that the agents are coupled according to the undirected cyclic 
graph shown in Fig. 3, since the force exerted by mass i on mass i + \ exactly reciprocates the force exerted 
by mass i -f 1 on mass i: 

ki+i,i{oti+i — oti) + ki+i,2{ai+i — tti) = —ki+i^i{ai — cti+i) — ki+i^2{<^i — oti+i)- 

We will treat the graph in Fig. 3 as a special case of directed graph with symmetric adjacency matrix. On 
the other hand, the communication topology of the system is assumed to be a directed graph shown in 
Fig 4. According to this graph, agents 1, 8, 12 and 15 observe the leader. 

To illustrate the results of the paper, the protocol matrices were computed using Theorem 2, and then 
the trajectories of the coupled mass-spring-damper system with the obtained protocol were simulated. To 
design the protocol, the parameters of the coupled mass-spring-damper system were chosen to be m = 3 
kg, k = 0.1 N/m, and k = 2.4 — lAp{t) N/m with p{t) = cos(f). The headway distance between the mass 
bodies was assumed to be 1 m. 

In the simulation, the spring constants and damper coefficients were chosen to be 1 =0.1 N/m and 
fci ,2 = 0.1 N/(m/s), 1 = 0,..., 20, respectively. The initial states of the leader and the followers were set 


© 

Figure 4: Directed communication graph. 
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to the values shown in table 1 . In the cost function, we chose Q 


and R = 0.001. 


10 0 
0 100 


Table 1: The initial states of the leader and followers. 


2;o(0) 

4(0) 

4(0) 

4(0) 

4(0) 

2^40) 

2^6(0) 

[0.5 0] 

[0.45 0] 

[0.4 0] 

[0.3 0] 

[0.2 0] 

[0.15 0] 

[0.25 0] 

4(0) 

4(0) 

4(0) 

2^10(0) 

xii(O) 

2^12(0) 

2^13(0) 

[0.35 0] 

[0.45 0] 

[0.55 0] 

[0.65 0] 

[0.55 0] 

[0.45 0] 

[0.35 0] 

44(0) 

2^15(0) 

440) 

2^17(0) 

2^18(0) 

2^19(0) 

2^20(0) 

[0.45 0] 

[0.5 0] 

[0.4 0] 

[0.3 0.1] 

[0.2 0.2] 

[0.1 0.3] 

[0 0.4] 


To design the synchronization protocol, we chose T = [—1,1] with 4 design points 

= { - 1, -0.3333,0.3333, l}. 

Choosing /Ig = 0.3111 so that the properties (37) hold for each design point, we solved the corresponding 
LMIs (18). Next, we constructed a continuous-gain control protocol by using the interpolation technique 
based on Theorem 2. We computed p = 0.0143, g = 0.0082 and we need to select q G (0,1) such 
that supt \p{t)\ < 1.7485(7, q G (0,1). Since supj |p(f)| = sup^ | — sinfj = 1, this allows us to choose 
q = 0.5750. Then the theoretically predicted bound on the performance (54) is computed to be 329.1316. 
On the other hand, by simulating the system on the time interval [0,12], we numerically found the value of 
the performance cost (6) for this particular instance of the system to be equal to 77 (m) = 39.7982. 

Moreover, switching control protocol which does not involve the interpolation technique was simulated 
in the example to compare difference between the gain scheduling and switching design. The computed 
performance cost (6) for the switching control of the system is equal to 77(w) = 40.2858. The simulation 
results are shown in Fig. 5, 6, 7, 8 and 9. 




Figure 5: Relative displacement (the top figure) and relative velocities of the followers with respect to the 
leader. 

The Fig. 5 demonstrates that the proposed continuous protocol enables all the followers to synchronize 
to the leader. The displacement and velocity trajectories of each block are shown in Fig. 6. Fig. 7 shows 
relative distances between the blocks ai — ai-i + h, where h — Im is the headway distance. From the 
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Figure 6; The displacement (the top figure) and velocity trajectories of the leader and followers. 



Figure 7: The relative distance between the blocks. 
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Figure 8: The continuous control signals (the top figure) and switching control signals of the followers. 
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Figure 9; Accelerations of the leader and the followers without control forces (the top figure), with contin¬ 
uous control forces and switching control forces (the bottom figure). 
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figure, the relative distances vary initially between 0.9m and 1.1m and remain positive. This indicates that 
collisions between the blocks were avoided. Fig. 8 confirms that the control signals are continuous for 
each follower by using the interpolation technique, while the switching control scheme leads to discontin¬ 
uous control signals. This also causes accelerations of the subsystems to be continuous and discontinuous, 
respectively; see Fig. 9. Discontinuous accelerations are often undesirable in practical systems. We also 
notice that the accelerations exerted by the proposed controllers are of comparable value with the acceler¬ 
ations occurring in the control-free system (Fig. 9); this indicates that the controller gains have been tuned 
to acceptable values. 

5 Conclusions 

The paper has considered the leader-follower control problem for a parameter varying system with 
directed communication topology and linear uncertain coupling, subject to norm-bounded constraints. In 
contrast to many existing works, we assume that the leader is selected among the agents, and remains 
coupled to the rest of the system. Therefore, in a sense the problem under consideration has addressed 
synchronization in leaderless interconnected multi-agent systems. The main technical challenge has been 
to overcome the effects of physical interconnections. 

For this problem, we have proposed a continuous gain-scheduled consensus-type control protocol 
which employs an interpolation technique. The sufficient condition for the synchronization of the sys¬ 
tem is obtained which guarantees a suboptimal bound on tracking performance. The condition is in the 
LMI form, and is numerically tractable, although it requires solving coupled LMIs. 
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